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Quantum Markov Process on a Lattice 
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We develop a systematic description of Weyl and Fano operators on a lattice phase space. Introducing the 
so-called ghost variable even on an odd lattice, odd and even lattices can be treated in a symmetric way. The 
Wigner function is defined using these operators on the quantum phase space, which can be interpreted as a spin 
phase space. If we extend the space with a dichotomic variable, a positive distribution function can be defined 
on the new space. It is shown that there exits a quantum Markov process on the extended space which describes 
the time evolution of the distribution function. 



1. Introduction 

The quantum Markov process for an integer 
spin system was constructed on the lattice phase 
space by Cohendet et al. PP. The time evo- 
lution of the pseudo-distribution function, i.e., 
the Wigner function on the space can be de- 
rived from the process. However, the extension 
of their method to a half integer spin system is 
not straightforward. There are difficulties in the 
construction of the Wigner function on an even 
lattice. A way to avoid this difficulties is to intro- 
duce the so-called " ghost variable" [2j . The Weyl 
and Fano operators can be constructed on an 
even lattice if we consider fictitious lattice points 
between the original ones (see also 3 ). In the 
present paper, we develop the systematic treat- 
ment of the Wigner function on both odd and 
even lattices and show that a quantum Markov 
process for a half integer spin system can also be 
constructed in a similar manner. 



2. Weyl and Fano operators with ghost 
variables on a lattice 

We consider a lattice composed of N lattice 
points. A lattice with odd or even lattice points 
is called odd or even lattice, respectively. We de- 
note the set of integers and half integers by Z. 
For odd N, we define J°,J° as the sets com- 
posed of integers, integers and half integers be- 
tween -(N/2) and (N/2) + l, respectively, and J£ 



be the non-negative part of J° . For even N, we 
define J e ,J e ,J+ similarly between —(N/2) and 
(N — l)/2. The symbols J, J, J+ are the abbrevi- 
ations of J°, J°, Jl for odd N and J e , J e , J% for 
even N . We consider a Kronecker's delta function 
Sn^ri on Z, defined by 5fn,l = (5 2m ,2n in mod 2N, 
where m, n G Z. We denote one of primitive roots 
of unity of order N by ui, e.g., w = exp (^p)- 

The (n, m) components of basic phase, shift 
and skew operators Q,P,T are defined as 

Q = (c-e„) ), P = (5%) lt J,T = (Cl,o), (1) 
for odd N, and 

Q = (w-jW), P = ((1 - 26^ )Ci J, 
T = (5 {N) ) 



(2) 



for even N, where m, n S J. The Froquet or 
Bloch angle is it in the latter case 0]. 

3. Wigner function on a lattice 

3.1. Weyl operators 

We define the Weyl operators by 



W„ 



-2mnn2n 



Q" l P 



-2m 



,2mn i 



P~ 2m Q Zn , (3) 



for both odd and even N, where m,n € J. We 
can see the Weyl operators satisfy the relations 
for reflection and shift: 



WL n = W_ m ,_ n = TW m ,„T, 



(4) 
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W W , / = cj- 2 ( mn '- nm ')w _l / 

" m,n yv m',n' — w vv m+m' 



n+n' • 



(5) 



The TV 2 operators W m ^ n with (to, n) in the 
basic region J + x J + are complete and orthog- 
onal in the trace norm. Generally, W m ,„ with 
(m,n) e Z x Z and W mod(m i f )imod( „ ,£) in the 
basic region are related as 
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(6) 



(8) 



for even N. 

The Weyl operators have the next symplec- 
tic property, i.e., they are rotated by ir/2 anti- 
clockwise under the Fourier transformation, 



W m . n = S(m,n)W mod{m K) tiaod ^ ni N), 
where 6(m,n) = (-l)^™."), and 

r(m, n) = 2n 

for odd N, 

r{m, n) = 2n 
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3.2. Fano operators 

We define the Fano operators by 



A m , n — W m . n T — ui 



-Iran Q2nrp p2m 



(9) 



(10) 



for both odd and even N, where to, n £ J. They 
are hermite but over complete. Those in the basic 
region are orthogonal and complete. The Fano 
operators have the same symplectic property as 
the Weyl operators, 



T A — A 



(11) 



We can see that the Fano operators have proper 
marginal properties: 

^ A m>n = 0, for half integer to (ghost) (12) 

ne.J 

A m ,„ = 2N$$5V2 ) , for integer to (13) 

ne.J 

^ A m>n = 0, for half integer n (ghost) (14) 

mGJ 



J2 A m,« = 2(u/™ ( ^' ) ), for integer n (15) 
for odd N , and 

^ A TOi „ = 0, for integer to (ghost) (16) 

neJ 



J2 A m ,„ = 2iV(<5^ ) <^ ) ),for half integer to(17) 

ne.J 

A TOi „ = 0, for integer n (ghost) (18) 

A mtn = 2(uj n ^), for half integer n (19) 
for even N. 

3.3. Time evolution of the Wigner function 

We define the Wigner function by 

W(m,n) = -^tr(pA ro ,„), (20) 
which is real valued and bounded by 1/N, 
-^<W(m,n)<± (21) 
The average value of an observable O is given by 
(0)=N 0(m,n)W(m,n), (22) 

where 0(m,n) = ^ tr(OA TOj „). We rewrite the 
time evolution equation of the density matrix p, 



rw l 



(23) 



in terms of the Wigner function. We expand the 
Hamiltonian iJ by using the Weyl operators in 
the basic region, and the density matrix by the 
Fano operators in the same region. Equating the 
coefficient of A m ,„, the time evolution equation 
of the Wigner function is given by 



dt 



-2 J2 H+(m",n") 



f 2ni 



sin j — {2(m"n' - ri'm!) - a(m" , n")} 
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(m" + m', n" + n')W m > n > 



(24) 



where m" = mod(m — m', y), i" = mod(n — 
n', y) and H + (m,n) and a(m, n) are deter- 
mined by the polar decomposition H(m, n) = 

4. Construction of a Markov process 

We extend the basic region with the dichotomic 
variable a € {±1} = B, i.e., J + x J + — > J + x 
J + x 5, and consider a new real valued function 
G(m, n, c) on the extended space, 

G(m,n,a) = ±(^ + aW(m,n)\. (25) 

It can be seen that G(to, n, cr) is positive and 
satisfies the inequality, 



f 3 

77 < G(m,n,cr) < 77, 

and normalization condition, 

^2 G(m,n,cr) = l. 

{m,n)e.J + X.J + , ctG{±1} 



(26) 



(27) 



The positivity of G(m, n, cr) follows from the 
boundedness of the Wigner function Ea. (|2T)l . The 
average value of an observable O is given by 

(0)=2N 2 (rO(m,n)G(m,n,a). (28) 

It is natural to call G(to, n, cr) as a distribution 
function on the extended space J+ x J + x B. The 
time evolution of G(m, n, a) is given by 



dt 



G(to, n, a) 



^2 -sgn(o-cr') 

m'n'eJ+,o-'e{±l} 



H + (m", n") sin I — {2(mV - n'W) - a(m", n")} 

x 0(m" + m', n" + n')G(m', ri , a'). (29) 
Introducing a generating operator as 

■A t (m, n, cr; m', n', cr') = H + (m", n") / — - 

|_G(m , n , cr) 



-sgn(crer') x sin ^ ^^{2(m"n' — n"m ) — a(m" , n")} 



6{m" + m',n" +n')}, 
for (m,n) ^ (m',n'), 
At{m, n, a; m, n, cr') = 0, 
for (m, n) = (m', n'), a 7^ cr', 

At{m, n, a; m, n, a) 



(30) 
(31) 



G(m, n, cr) 



5^ W + (m',n') 

(m',n')eJ+xJ+\(0,0) 



(32) 



for (m,n,a) = (ml ', n' , cr') , the equation for 
G(to, n, cr) can be written briefly as 



(ft 



G(m, n, cr) 



E 



rn>,n'£j + , ctG{±1} 

At(m, n, cr; m', n', (j')G(m' ', n', cr'). 



(33) 



It can be checked that the generator 
At(m,n,a;m',n',a') satisfies the Markov con- 
dition, 

V] A(m,n,a;m',n',a') = 0, (34) 

(m,n,(r)£j + X J + xB 

^4.(m, n, cr; m , n , cr') > 0, if (to, n, cr) ^ (to', n, cr'), (35) 

which assures the existence of a background quan- 
tum Markov process on the extended lattice 
quantum phase space. 

REFERENCES 

1. O. Cohendet, Ph. Combe, M. Sirugue, 
M. Sirugue-Collin, J. Phys. A : Math. Gen. 
21 (1988) 2875-2883. 

2. U. Leonhardt, Phys. Rev. Lett. 74 (1995) 
4101; Phys. Rev. A 53 (1996) 2998. 

3. A. Takami, T. Hashimoto, M. Horibe, 
A. Hayashi, Phys. Rev. 64 (2001) 032114, 1- 
6; M. Horibe, A. Takami, T. Hashimoto, 
A. Hayashi, Phys. Rev. 65 (2002) 032105, 1-6. 

4. A. M. F. Rivas, A. M. Ozorio de Almeida, 
Annals Phys. 276 (1999) 223-256. 



